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Abstract 

Background: Murine leukemia virus-type endogenous retroviruses (MuLV-ERVs) constitute -10% of the mouse 
genome and are associated with various pathophysiologic processes. In this study, we examined whether 
MuLV-ERVs' response to burn-elicited stressors is specific for certain lymphocyte populations and/or locations of 
lymphoid organ. 

Results: B- and T-cells, which were sorted from nine lymphoid organs of C57BL/6J mice after burn, were subjected 
to MuLV-ERV expression analyses. Overall, the post-burn MuLV-ERV expression pattern was dependent on 
lymphocyte type, time after injury, location of lymphoid organ, and MuLV-ERV type. For instance, the MuLV-ERV 
expression in T-cells from the thymus and three cervical lymph nodes decreased at 3 hours post-burn while the 
expression of some MuLV-ERVs was augmented in B-cells derived from the mesenteric lymph node. The MuLV-ERV 
U3 sequences population of the burn-24 hours group was less diverse in comparison to the no burn and 
burn-3 hours groups. In addition, it was apparent that at the 24 hours time point, the U3 populations of B-cells 
from both no burn and burn groups were less heterogeneous than the T-cells' U3 populations. Using the U3 
sequences, some of which were isolated only from specific experimental groups (B- vs. T-cells; no burn vs. burn), as 
probes, 51 putative MuLV-ERVs, including 16 full-length proviruses, were mapped followed by characterization of 
their biologic properties. 

Conclusion: MuLV-ERVs' response to burn-elicited stressors may be differentially controlled depending on 
lymphocyte type, location of lymphoid organ, MuLV-ERV type, and stress duration. 
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Background 

Proviral sequences of the murine leukemia virus-type 
endogenous retroviruses (MuLV-ERVs) are genomic 
footprints reflecting the ancient infection of germ-line 
cells and are presumed to constitute ~10% of the mouse 
genome [1]. Although most MuLV-ERV genomes are 
reported to be defective primarily due to mutation and 
recombination events, some retain intact proviral gen- 
ome features, such as identical long terminal repeats 
(LTRs) at both ends, which are formed by retroviruses' 
unique replication process, and coding potential for the 
major structural genes (gag, pol and env) which are 
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essential for virion assembly [2-4]. Expression of MuLV- 
ERVs is mostly controlled by their U3 promoters located 
on the 5' LTR, in conjunction with the transcription ma- 
chinery that is dynamically and uniquely formed within 
individual cells [5]. 

Recent studies provide evidence that human endogen- 
ous retroviruses (HERVs) are involved in normal physi- 
ology as well as various disease processes [6]. The 
HERV-W envelope protein, called syncytin, is highly 
expressed in the syncytiotrophoblast layer and plays a 
critical role in placental differentiation [7]. In addition, 
the inflammatory properties of syncytin contribute to 
the degeneration of oligodendrocytes responsible for 
the development of multiple sclerosis, an autoimmune 
disease [8]. 
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Patients with severe burn injury often succumb to sys- 
temic immune disorder in conjunction with internal 
organ damage [9]. However, the mechanisms underlying 
the post-burn pathologic changes have not yet been fully 
characterized. The results from recent studies suggest 
that post-burn systemic inflammatory response is asso- 
ciated with the activation and/or depletion of specific 
lymphocyte subpopulations, such as CD4 + CD25 + T regu- 
latory cells, lymphocytes expressing CD25, CD69 and 
CD71, NK-T-cells, and CD8 + T-cells [10-12]. Sepsis- 
elicited stressors induce apoptosis in the spleen, resulting 
in the depletion of B-cells and accompanying immune 
malfunction [13]. However, the central signaling mole- 
cules responsible for the post-burn activation and/or de- 
pletion of lymphocytes remain to be elucidated. 

Our recent studies demonstrated that burn-elicited 
stressors differentially alter the expression of MuLV-ERVs 
in a range of internal organs, including various lymphoid 
tissues, in a ERV type- and organ-specific manner [14]. 
Moreover, some MuLV-ERV envelope proteins were able 
to modulate the mRNA expression of inflammatory cyto- 
kines such as IL-6 and IL-lp in macrophage cells [15]. In 
this study, we hypothesize that the differential post-burn 
alterations in MuLV-ERV expression, depending on 
lymphocyte types and/or lymphoid organs, may be asso- 
ciated with systemic pathogeneses, including immune dis- 
order. We examined whether MuLV-ERVs' response to 
burn-elicited stressors is specific depending on lympho- 
cyte type and the location of lymphoid organ followed by 
an investigation of the biologic characteristics of the rele- 
vant MuLV-ERVs. 

Methods 

Animal experiments 

The experimental protocols involving mice were approved 
by the Animal Use and Care Administrative Advisory 
Committee of the University of California, Davis. Female 
C57BL/6J mice (12 weeks old) were purchased from 
Jackson Laboratory- West (Sacramento, CA). The baseline 
MuLV-ERV expression patterns were examined in nine 
lymphoid organs (spleen, thymus, and seven lymph nodes 
[axillary, inguinal, mesenteric and pooled cervicals {man- 
dibular, accessory mandibular, superficial parotid, and 
deep cervical}]). Three mice were sacrificed for tissue col- 
lection by cervical dislocation without any treatment. The 
burn protocol has been described previously [16]. Briefly, 
under general anesthesia, an ~18% total body surface area 
flame burn was generated on the shaved back of mice fol- 
lowed by resuscitation with saline and buprenorphine. 
Control (no burn) mice were treated the same as burn 
mice but without the flame burn. The three mice from 
each experimental group were sacrificed by CO2 inhal- 
ation at 3 hours and 24 hours after burn followed by the 
collection of the same nine lymphoid organs as above 



without any pooling. The burn experiment was repeated 
four times. 

Cell sorting 

The collected tissues were minced, suspended in a tryp- 
sin/EDTA solution, and then passed through a cell 
strainer (BD Biosciences, CA). The cells were washed 
with PBS containing 0.1% BSA followed by centrifuga- 
tion at 1,000 xg for 5 min at 4°C. Following the removal 
of red blood cells, cells were sorted into B- and T-cells 
using pan-B and pan-T Dynabeads (Invitrogen, Carlsbad, 
CA), respectively, according to the protocols provided 
by the manufacturer. 

Total RNA isolation and RT-PCR 

Total RNAs were isolated from the sorted cells using an 
RNeasy Mini Kit (Qiagen, Valencia, CA) following 
homogenization using a QIAshredder (Qiagen). cDNAs 
were synthesized using Sensiscript reverse transcriptase 
(Qiagen), RNase inhibitor (Promega, Madison, WI) and 
an oligo-dT primer (5'-GGC CAC GCG TCG ACT 
AGT ACT TTT TTT TTT TTT TTT T- 3'). A set of 
primers, ERV-U1 (5'-CGG GCG ACT CAG TCT ATC 
GG-3') and ERV-U2 (5'-CAG TAT CAC CAA CTC 
AAA TC-3'), was used to amplify the MuLV-ERV U3 
regions [17]. P-actin, serving as a normalization control, 
was amplified with the primer set: 5'-CCA ACT GGG 
ACG ACA TGG AA-3' and 5'-GTA GAT GGG CAC 
AGT GTG GG-3'. 

Measurement of relative purity of sorted cell populations 

The relative purity of the sorted cell populations was 
measured by real-time RT-PCR using CD20 (B-cell anti- 
gen) primers (5'- AAA ACC TCC AGG AAG AGT TTG 
GTC-3' and 5'-CGA TCT CAT TTT CCA CTG GCA 
AG-3') for B-cells and Thyl (Thymus cell antigen-1) pri- 
mers (5'-TCC TGC TCT CAG TCT TGC-3' and 5'- 
AGT AGT CGC CCT CAT CCT TG-3') for T-cells. 
Real-time RT-PCR was performed using the Brilliant-III 
SYBR Green QPCR kit from Agilent Technologies (La 
Jolla, CA) and the results were analyzed using the 
MxPro (v.4.10) program (Agilent Technologies). Cell 
type specific mRNA levels (CD20 and Thyl) were calcu- 
lated from each sorted cell population with the level of 
[3-actin as a normalization control, using the same pri- 
mer set as above. The relative purity of each cell popula- 
tion was calculated with the assumption that there are 
only B- and T-cells in each sorted cell population and 
CD20 and Thyl are exclusively expressed in B-cells and 
T-cells, respectively. 

Densitometry and statistical analysis 

The density of the PCR products was measured using 
KODAK Molecular Imaging Software Version 4.5.1 
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(Eastman Kodak Company, Rochester, NY). All optical 
densities were measured based on mass values with the 
Hi-Lo DNA marker (a gift from Minnesota Molecular 
Inc., Minneapolis, MN). The density of each MuLV-ERV 
U3 band was normalized to its corresponding (3-actin 
value to obtain the relative intensity of each U3 band. 
Differences in the relative intensity of the matching 
MuLV-ERV U3 bands between no burn and burn were 
calculated as fold changes. Statistical significance (p < 
0.05) of fold changes between the matching pairs (no 
burn and burn) of the individual U3 bands was deter- 
mined using the data from four burn experiments and 
Student's i-test. 

Cloning and sequencing 

Eighty-eight MuLV-ERV U3 PCR amplicons/bands, 
which had contrasting densities between individual no 
burn and burn pairs, were selected for sequence ana- 
lyses. The PCR amplicons were purified using the QIA- 
quick Gel Extraction kit (Qiagen) and then cloned into 
the pGEM®-T Easy vector (Promega). Three clones 
were prepared from each U3 PCR amplicon/band. The 
plasmid DNAs were prepared using a QIAprep Spin 
Miniprep kit (Qiagen) for sequencing analysis. Sequen- 
cing was performed at a core sequencing facility at the 
University of California, Davis. 

Alignment and phylogenetic analyses 

From a total of 260 U3 sequences, 100 unique U3 
sequences were identified using Seqman Pro Version 8.1.2 
(DNASTAR, Madison, WI). Multiple sequence alignment 
of the 100 unique U3s was performed using the Clustal W 
method in MegAlign (DNASTAR), and a phylogenetic 
tree was generated by the HYPERTREE program [18]. 
Seven U3 sequences, representing the seven phylogenetic 
branches, were subjected to a multiple alignment analysis 
to examine the U3 sequence features of direct repeats (1/ 
1*, 4/4*, 5/5*, and 6/6*), insertion of 190 nucleotides, and 
unique region (2) based on the protocol by Tomonaga 
etal. [17,19]. 

Profiling of transcription regulatory elements 

The 100 unique MuLV-ERV U3 sequences were surveyed 
for putative transcription regulatory elements using the 
Matlnspector program (Genomatix, Munich, Germany) 
within the vertebrate matrix group with a core similarity 
higher than 0.95 [20]. 

Determination of population diversity of U3 sequences 
within each experimental group 

Using the MEGA4 program [21], the population diver- 
sity of U3 sequences was measured in each experimental 
group and compared among different groups in two ana- 
lyses: a) no burn, burn-3 hours, and burn-24 hours, and 



b) no burn-B-cells, burn-24 hours-B-cells, no burn-T- 
cells, and burn-24 hours-T-cells. 

In silico mapping of putative MuLV-ERVs and 
characterization of their biological properties 

The C57BL/6J mouse genome database from the Na- 
tional Center for Biotechnology Information (NCBI) was 
surveyed for putative MuLV-ERVs with a greater than 
98% identity using the 100 unique U3 sequences as min- 
ing probes and the Advanced BLAST program (NCBI). 
The coding potential of each MuLV-ERV was deter- 
mined using the Vector NTI program (Invitrogen). 
Three polypeptide genes {gag, pol, and env) were noted 
as intact (+), partial (P), or defective (-) relative to the 
references (AY219567.2 and AF033811). Each MuLV- 
ERV was mapped with the chromosomal location and 
strand orientation. The primer binding site (PBS) of each 
MuLV-ERV was identified by examination of a stretch of 
18 nucleotides immediately downstream of the 5' LTR 
region [22]. In addition, the annotated gene(s) within 
100 Kb upstream and downstream of each MuLV-ERV 
locus, if any, were mapped. 

Results and discussion 

Lymphocyte type-specific MuLV-ERV expression patterns 
in lymphoid organs of normal mice 

Prior to the investigation of post-burn alterations in 
MuLV-ERV expression in B- and T-cells of different 
lymphoid organs (spleen, thymus, axillary lymph node 
[LN], inguinal LN, mesenteric LN, and cervical LNs), 
baseline MuLV-ERV expression levels were examined in 
B- and T-cells isolated from -12 week old female 
C57BL/6J mice. The relative purity of the sorted cell 
populations was determined by measuring CD20 (for B- 
cells) and Thyl (for T-cells) levels by real-time RT-PCR. 
The relative purity of the sorted B- and T-cell popula- 
tions was greater than 93%, except for the B-cell popu- 
lation in the thymus (-65%). Throughout the lymphoid 
organs examined, the MuLV-ERV expression pattern, 
which was examined by RT-PCR amplification of the 3' 
U3 regions, was substantially different between B- and 
T-cells (Figure 1). The expression pattern of an intense 
-500 bp and two (-750 bp and -450 bp) lighter U3 
bands was shared by B-cells derived from all lymphoid 
organs examined. In contrast, four distinct U3 bands, 
ranging from -500 bp to -750 bp, were amplified in T- 
cells with variable intensity depending on the source 
lymphoid organs. It was interesting to note that the 
-500 bp U3 band was more intense in B-cells from all 
four LN samples compared to B-cells from the spleen 
and thymus. The findings from this study suggest that 
the expression of MuLV-ERVs is differentially regulated 
depending on the lymphocyte type, source lymphoid 
organ, and ERV type. 



Lee ef al. BMC Immunology 2013, 14:2 
http://www.biomedcentral.eom/1471-2172/14/2 



Page 4 of 10 



pol 



3'LTR ( 



Cell type 



E 5, 

> 0) 



ERVU3 



B-actln 



B 



SM OT 



= C 



■= S » 



x z °>z » ; 
<-i=_i 5 




SM 



a 



= c 
E ™ 

>» ^ 



<fl CO) 

o z o z 
E _i O _i 




750 



-500 

_300 

-200 
(bp) 



Figure 1 MuLV-ERV expression in B- and T-cells derived from lymphoid organs of C57BL/6J mice. MuLV-ERV expression was examined in 
B- and T-cells isolated from the spleen, thymus, and seven different lymph nodes of female C57BL/6J mice without any treatment. A schematic 
drawing of a typical MuLV-ERV proviral structure indicates the relative locations of a set of primers used for amplification of the U3 regions. 
Differential MuLV-ERV expression patterns were observed in B- and T-cells derived from various lymphoid organs. LN (lymph node), SM (size marker). 



Burn-elicited alterations in MuLV-ERV expression in 
lymphocytes depending on lymphocyte type, source 
lymphoid organ, and ERV type 

The effects of burn-elicited stressors on the expression of 
MuLV-ERVs in the B- and T-cells from the spleen, thymus, 
and seven lymph nodes were examined at 3 hours and 
24 hours after injury by amplifying the 3' U3 regions. At 
3 hours post-burn, the MuLV-ERV expression, represented 
by four distinct U3 bands, was decreased in the T-cells from 



the thymus and three cervical LNs (mandibular, accessory 
mandibular, and superficial parotid), but not in the T-cells 
from the deep cervical LN (Figure 2). In contrast, there was 
an increase in intensity of one to two U3 bands in the T- 
cells from the thymus and superficial parotid cervical LN at 
24 hours post-burn when there was a decrease in U3 bands 
in the T-cells of the spleen and mesenteric LN. 

In B-cells, there was an induction of two (a and b 
bands) of the three U3 bands in the mesenteric LN at 
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Figure 2 Post-burn alterations in the expression of MuLV-ERVs in B- and T-cells isolated from lymphoid organs of C57BL/6J mice. Post- 
burn (3 hours and 24 hours) changes in the expression of MuLV-ERVs were measured in the B- and T-cells isolated from nine different lymphoid 
organs. The results, which were presented as a fold change (burn/no burn) for the individual MuLV-ERV U3 bands (a, b, and c for B-cells; d, e, f, 
and g for T-cells), revealed a unique pattern for each experimental group (cell/organ/time). SPL (spleen), THY (thymus), LN (lymph node), AXI 
(axillary LN), ING (inguinal LN), MES (mesenteric LN), MAN (mandibular cervical LN), ACC (accessory mandibular cervical LN), SUP (superficial 
parotid cervical LN), DEE (deep cervical LN). A black arrow indicates unmeasured samples due to the low expression level. Error bars indicate 
standard error of the mean. * P < 0.05. 
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3 hours post-burn while two bands (a and c) were down 
regulated in the B-cells of the mandibular cervical LN. 
However, at 24 hours post-burn, the expression levels 
of all three U3 bands were increased in the B-cells from 
the mandibular cervical LN. Also at the 24 hours time 
point, post-burn increases in one or two U3 bands were 
also observed in the B-cells of the spleen, thymus, axil- 
lary LN, inguinal LN, and/or accessory mandibular 
cervical LN. 

The results from this study suggest that the lymphocyte 
type (B- or T-cell), source lymphoid organ, and duration 
after injury play a role in post-burn changes (induction or 
repression) in the expression of MuLV-ERVs. It is possible 
that the post-burn differential expression of MuLV-ERVs 
resulted from changes in the composition of the B- and 
T-cell sub-populations (e.g. naive vs. effector cells), which 
have different transcription environments, instead of 
(and/or in addition to) the affected MuLV-ERVs' altered 
transcriptional potentials. It was interesting to note that 
the post-burn increase in the MuLV-ERV expression in 
the B-cells from the mesenteric LN parallels the prolifera- 
tion of specific B-cell populations in the mesenteric LN 
and spleen after burn [23]. 

MuLV-ERV U3 sequences identified in specific lymphocyte 
type and/or injury (no burn/burn) condition 

Three clones from each of the 88 MuLV-ERV U3 ampli- 
cons/bands, which showed contrasting expression levels 
between the individual no burn and burn pairs, were 
sequenced to identify lymphocyte type- and/or injury 
(no burn/burn)-specific U3 sequences. From a total of 
260 MuLV-ERV U3 sequences cloned from the entire set 
of MuLV-ERV U3 bands examined, 100 unique U3 
sequences, ranging from 346 nucleotides to 615 nucleo- 
tides in size, were identified (Additional file 1: Table SI). 
An initial phylogenetic analysis resulted in seven distinct 
branches of U3 sequences, primarily segregated by their 
size (group I [largest U3s] ~ group VII [smallest U3s]) 
(Figure 3). Group III (440 bp length) had the most num- 
ber of U3 sequences (40) followed by group VII (346 bp 
length) (25 sequences) then group I (greater than 603 bp 
length) (19 sequences). There were only two U3 sequences 
in groups IV and VI (433 bp and 361 bp in size, respect- 
ively). The group II U3 sequences (5) were isolated only 
from T-cells and the group VI U3 sequences (2) were 
unique for B-cells. On the other hand, the group IV U3 
sequences (2) were identified only from no burn control 
mice. The other four branches (groups I, III, V, and VII) 
consisted of U3 sequences isolated from both lymphocyte 
types derived from both no burn and burn mice. 

The results obtained from this study suggest that some 
MuLV-ERVs are selectively expressed depending on 
lymphocyte type (B- or T-cell) and/or presence of burn- 
elicited stressors. However, it is possible that an expanded 
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Figure 3 Phylogenetic relatedness of 100 unique MuLV-ERV U3 
sequences. Phylogenetic relatedness of the 100 unique MuLV-ERV 
U3 sequences was determined by a multiple alignment analysis. The 
MuLV-ERV 113 sequences in each branch were analyzed for their 
presence in a specific lymphocyte type (B- or T-cell) and/or stress 
condition (no burn or burn). Seven branches were formed, primarily 
based on the size of the U3 sequences. Number of sequences is 
in parentheses. 



saturation cloning protocol would find that some cell type 
and/or burn-specific U3 sequences are not exclusively 
present in one group over the other. In this case, a statis- 
tical evaluation may be needed to determine the specificity 
of the U3 sequences in regard to the lymphocyte type 
and/or burn-elicited stressors. 

Characteristics of seven U3 sequences representing the 
phylogenetic branches 

Seven representative U3 sequences, which were chosen 
from the seven main phylogenetic branches, were subjected 
to an alignment analysis using reference sequences to 
characterize their sequence features, such as direct repeats, 
unique sequence, and an insertion of 190 nucleotides. 
These U3 sequence features were utilized as key parameters 
for the in silico determination of the MuLV-ERV tropism 
traits [19]. The U3 sequences of the groups I and II had an 
insertion of 190 nucleotides, which is consistent with their 
relatively large size compared to the other five branches 
(Table 1). There was a well-conserved TATA box in all 
seven U3 sequences; however, each U3 sequence retained a 
unique profile in regard to the direct repeats (1/1*, 4/4*, 
5/5*, and 6/6*) and the unique sequence (2). The specific 
features of the individual U3 sequences were translated into 
the tropism traits of the respective putative MuLV-ERVs 
[17,19]. However, the sequences features of the group VI 
U3 sequences were too diverse to predict its tropism trait. 
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Table 1 Sequence characteristics and putative tropism traits of seven MuLV-ERV U3 sequences representing the seven 
main phylogenetic branches 



Group 



U3 Probe 



Size (bp) 



1/1* 



Repeat/unique region 

2 4/4* 5/5* 



6/6* 



TATA 



Insertion 
(190 bp) 



Tropism 



B3h-Tth _ 1 

B3H"T|N"3 



603 
556 



+/+ 
+/+ 



+/+ 
+/+ 



D 3H~ ' TH~ 



440 



+/+ 



+ 



+/+ 



+/+ 



X-ll 



IV 
V 



NB 24H -T sc -3 
B3H-TTH-6 



433 
392 



+/+ 
+/+ 



+/+ 
+ 



+/+ 
+/+ 



+ 

+/+ 



X-ll 
X-l 



361 



+/+ 



+/+ 



+/+ 



VII 



i-Tth-1 1 



346 



+/+ 



+/+ 



+/+ 



Sequence characteristics: direct repeats (1/1*, 4/4*, 5/5*, and 6/6*), unique region (2), insertion of 190 nucleotides, and TATA box; +/+ (direct repeat), + (single 
copy), • (absence). 



Evidence for post-burn increase in the expression of 
specific MuLV-ERVs 

The population diversity of the U3 sequences (a total of 
260) was measured in each experimental group and 
compared among the groups in two different analyses 
[21]: i) no burn, burn-3 hours, and burn-24 hours, and 
ii) no burn-B-cells, burn-24 hours-B-cells, no burn-T- 
cells, and burn-24 hours-T-cells. The results from the 
first diversity analysis revealed that the U3 population of 
the burn-24 hours group was substantially less heteroge- 
neous compared to the other two groups (no burn and 
burn-3 hours) (Figure 4). In contrast, it appears that 
there was an increase in diversity among the U3 popula- 
tion in the burn-3 hours group. The finding of decreased 
U3 population diversity within the burn-24 hours group 
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Figure 4 Decreased population diversity of MuLV-ERV U3 
sequences at 24 hours after burn. The U3 population diversity 
within each experimental group was measured [21] and the relevant 
values, indicating different diversity levels, were presented as 
nucleotide diversity and Tajima statistics scores. There was a marked 
decrease in the U3 population diversity at 24 hours after burn 
compared to no burn. 
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Table 2 Differential transcription potentials of MuLV-ERV 
U3-promoter sequences 

Occurrence 
Burn/No burn 



DMRT 


ctttcttGTAActtttgagaa 


68 


/ 


1 


BRNF 


tctggatctatTGATttga 


35 


/ 


A 


BCL6 


ctctttcTGGAaactga 


31 


/ 


3 


TF2D 


gctctggtactttttcatcccttgcaaaatGGCGttact 


27 


/ 


1 


PTF1 


gccaTCTGttcttggccctga 


28 


/ 


1 


MYT1 


tgagAGCTcagct 


27 


/ 


2 


AIRE 


gtatttctcggtcatttGGGGaaactg 


19 


/ 


0 


ZNFP 


cacaCACCtggtc 


22 


/ 


0 


NRF1 


gggGCCCaggcgcttga 


21 


/ 


1 


TRE 


Consensus motif 


Occurrence 


B-cell/T-cell 


PAX2 


gttattcggggaacctgAAACtg 


32 


/ 


13 


HAML 


ctgaGTGGttagttc 


35 


/ 


8 


NKXH 


acttttgaGAACttagctc 


33 


/ 


10 


HOMF 


taatctTTAActgttaatc 


16 


/ 


3/ 


SORY 


gttctaaaCAATgtctattttca 


14 


/ 


28 


HEAT 


gagagggtctaatGAAGgttcaaca 


10 


/ 


30 


MAZF 


atggGAGGggtac/gtacGAGGagagg 


8 


/ 


18 


MEF2 


acaatgtCTATtttcaagaaatg 


5 


/ 


19 


GATA 


tttAGATctaaac 


6 


/ 


16 


CART 


gggtcTAATgaaggttcaaca 


8 


/ 


16 


HOXF 


gggtcTAATgaaggttcaa 


7 


/ 


1/ 


LHXF 


gagagggtcTAATgaaggttcaa 


7 


/ 


16 


FXRE 


AGGTgcttgacca 


5 


/ 


19 


GRHL 


cacacaGGTTcaa 


7 


/ 


16 


NKX6 


atctTTAActgttaa 


6 


/ 


16 


CEBP 


aaggttcaGCAAtgg 


6 


/ 


lb 


FKHD 


tagttcTAAAcaatgtc 


5 


/ 


16 


PDX1 


gaaagagtcTAATgaaggt 


5 


/ 


15 


RXRF 


cacagGTTCaaggagtggccagagc 


5 


/ 


16 



Numbers indicate the occurrence frequency of specific TREs (transcription 
regulatory elements) within the 100 unique 113-promoter sequences. Please 
refer to (Additional file 3: Table S3) for the complete list of TRE abbreviations. 
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Table 3 Genomic map and characteristics of putative MuLV-ERVs harboring the unique U3 sequences 



No 


Virus 


Chr 


Location 


Ori 


Virus 


PBS 




ORF 




Occurrence 


Reference 


From 


To 


size (bp) 


gag 


pol 


env 


frequency 


1 


MuERV-49 


1 


133539144 


133548197 


+ 


9054 


Q 


- 


+ 


+ 


UAJ 


1,11,111,1V 


2 


MuERV-31 


1 


166153938 


166159683 


+ 


5746 


Q 


+ 


- 


- 


U/T 


1,11,1V 


3 


MuERV-30 


1 


184187251 


184196231 


+ 


8981 


Q 


+ 


P 


+ 


U/T 


1,11,1V 


4 


MuERV-33 


1 


193757492 


193766472 


+ 


8981 


Q 


+ 


P 


- 


U/T 


1,11,1V 


5 


MuERV-35 


2 


57066823 


57075803 


- 


8981 


Q 


+ 


+ 


+ 


U/T 


1,11,1V 


6 


MuERV-3 


2 


156182292 


156189341 


- 


7050 


Q 


+ 


- 


+ 


U/U 


II l,IV 


7 


MuERV-48 


3 


66900000 


66907716 


+ 


7717 


Q 


+ 


- 


- 


U/U 


1,11,111,1V 


8 


MuERV-55 


3 


151986104 


151995146 


+ 


9043 


Q 


+ 


+ 


+ 


U/U 


1,11,111,1V 


9 


MuERV-29 


4 


101535762 


101544743 


- 


8982 


Q 


+ 


P 


+ 


U/T 


1,11,1V 


10 


MuERV-28 


4 


107818639 


107827621 


- 


8983 


Q 


+ 


+ 


+ 


U/T 


1,11,1V 


11 


MuERV-2 


4 


132651872 


132657538 


+ 


5667 


Q 


- 


- 


- 


U/U 


l,ll,lll,IV 


12 


MuERV-6 


4 


133715306 


133720617 


+ 


5312 


Q 


+ 


- 


- 


U/U 


l,ll,lll,IV 


13 


MuLV_ERV 4 _ 5 * 


4 


145686198 


145694956 


- 


8759 


Q 


+ 


+ 


+ 


U/U 


this study 


14 


XMV9 


4 


146592937 


146601658 


+ 


8722 


Q 


+ 


+ 


+ 


U/U 


V 


15 


L-3-9.5 


5 


23206397 


23215063 


- 


8667 


P 


+ 


+ 


+ 


U/U 


1,11 


16 


MuLV_ERV 5 - 2 * 


5 


23721 126 


23727971 


- 


6846 


Q 


+ 


- 


+ 


U/U 


this study 


17 


MuERV-27 


5 


33864636 


33873616 


+ 


8981 


Q 


- 


+ 


+ 


U/T 


1,11,1V 


18 


MuERV-18 


5 


44525599 


44534580 


+ 


8982 


Q 


- 


+ 


+ 


U/T 


1,11,1V 


19 


MuERV-19 


5 


77359541 


77368522 


+ 


8982 


Q 


+ 


+ 


+ 


U/T 


1,11,1V 


20 


MuERV-26 


7 


7022622 


7031603 


- 


8982 


Q 


P 


- 


- 


U/T 


1,11,1V 


21 


MuERV-24 


7 


30400226 


30408924 


+ 


8699 


Q 


+ 


+ 


P 


U/T 


1,11,1V 


22 


MuERV-25 


7 


31473659 


31482639 


+ 


8981 


Q 


+ 


P 


- 


U/T 


1,11,1V 


23 


MuERV-34 


7 


123575623 


123581296 


+ 


5674 


Q 


+ 


- 


- 


U/T 


1,11,1V 


24 


MuERV-10 


8 


44397084 


44405817 


- 


8734 


Q 


+ 


+ 


+ 


U/T 


1,11,1V 


25 


MuERV-1 1 


8 


87649754 


87656141 


- 


6388 


Q 


+ 


P 


- 


U/T 


1,11,1V 


26 


MuERV-1 


8 


93410190 


93415857 


+ 


5668 


Q 


- 


- 


- 


U/U 


l,ll,lll,IV 


27 


MuERV-12 


8 


122413465 


1 22420940 


- 


7476 


Q 


- 


- 


+ 


U/T 


1,11,1V 


28 


MuERV-8 


9 


41713568 


41720351 


- 


6784 


Q 


- 


+ 


+ 


U/U 


1,11,1V 


29 


iLN(lll)-1a 


9 


62286981 


62295739 


+ 


8759 


Q 


+ 


+ 


+ 


U/U 


1 


30 


K-4-1 1.10a 


10 


8263933 


8271066 


- 


7134 


Q 


- 


- 


- 


U/T 


1 II 


31 


MuERV-53 


10 


22454084 


22463128 


+ 


9045 


Q 


+ 


+ 


+ 


U/U 


1,11,111,1V 


32 


MuERV-17 


10 


41 156454 


41 165434 


+ 


8981 


Q 


+ 


+ 


+ 


U/T 


1,11,1V 


33 


BM-a-Zlle 


11 


5905655 


591 1518 


+ 


5864 


Q 


- 


- 


- 


B/B 


III 


34 


MuERV-22 


11 


6649425 


6658394 


- 


8970 


Q 


+ 


- 


+ 


U/T 


1,11,1V 


35 


MuERV-23 


11 


820301 


8829280 


+ 


8980 


Q 


+ 


+ 


+ 


U/T 


1,11,1V 


36 


MuERV-50 


11 


60392151 


603995 1 0 


- 


7360 


Q 


- 


P 


- 


U/U 


1,11,111,1V 


37 


MuERV-51 


11 


76362310 


76371354 


+ 


9045 


Q 


+ 


+ 


+ 


U/U 


l,ll,lll,IV 


38 


MuERV-38 


11 


88751494 


88758280 


- 


6787 


Q 


+ 


- 


+ 


U/T 


1,11,1V 


39 


MuERV-37 


12 


20704753 


20710513 


+ 


5761 


Q 








U/T 


1,11,1V 


40 


MuERV-16 


12 


55875211 


55884192 


+ 


8982 


Q 


+ 


+ 


+ 


U/T 


IV 


41 


MuERV-5 


13 


67985825 


67994511 




8687 


Q 


+ 


+ 


+ 


U/U 


1,11,111,1V 


42 


MuERV-21 


13 


99756230 


99765210 




8981 


Q 


+ 


P 


+ 


U/T 


1,11,1V 


43 


MuERV-7 


14 


5491 2048 


54917591 




5544 


Q 


+ 






U/U 


IV 


44 


L-1-2.14 


14 


55149203 


55157025 




7823 


Q 


P 


+ 


+ 


U/U 


1,11 


45 


MuERV-15 


15 


76390247 


76397432 




7186 


Q 






+ 


U/T 


1,11,1V 


46 


U-1-5.16 


16 


76135288 


76144269 




8982 


Q 




P 




U/T 


1,11 
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Table 3 Genomic map and characteristics of putative MuLV-ERVs harboring the unique U3 sequences (Continued) 

47 MuERV-14 16 93697702 93706682 8981 Q + + U/T l,IV 

48 MuERV-20 18 82862943 82871822 + 8880 Q + + U/T 1,11,1V 

49 MuERV-32 19 38449512 38456573 + 7062 Q P U/T IJI.IV 

50 MuERV-9 19 61001210 61009936 - 8727 Q + + + U/U l,ll,IV 

51 MuERV-13 X 15051174 15060154 + 8981 Q + P + U/T l,ll,IV 

Chr (chromosome), Ori (orientation), PBS (primer binding site), Q (tRNA 9lutamir " ! ), P (tRNA proline ), ORF (open reading frame), + (intact/full-length), - (defective), U 
(both B- and T- cells; both no burn and burn), B (B-cell), T (T-cell), references I [14], II [33], HI [31], IV [30] and V [32]. *MuLV-ERV newly identified in this study. The 
names of MuLV-ERVs, which have been reported in previous studies, are adapted from the relevant reference in bold. 



was confirmed by its negative score for the Tajima's Neu- 
trality test in comparison to the positive values of the 
other two groups [24]. The results from the second di- 
versity analysis demonstrated that there is a decrease in 
U3 population diversity in both burn-24 hours-B-cells 
and burn-24 hours-T-cells groups in comparison to the 
no burn controls. Furthermore, it was apparent that at 
the 24 hours time point, the U3 population in B-cells 
from both no burn and burn were less diverse in com- 
parison to T-cells' U3 population. 

The findings from these diversity analyses indicate that 
burn-elicited stressors lead to the induction of specific 
MuLV-ERVs in a time after injury-dependent manner, 
and the population of MuLV-ERVs expressed in T-cells 
is more heterogeneous than the population in B-cells. 

Differential transcription potentials of MuLV-ERV 
U3-promoter sequences 

The U3 sequence on the 5' LTR, which is presumed to 
be identical to the U3 of the 3 LTR, serves as an ERV 
promoter. In this study, the transcription potentials of 
the 100 unique U3-promoter sequences were surveyed 
for transcription regulatory elements (TREs). Nine TREs 
were predominantly present in the U3 sequences isolated 
from both B- and T-cells of burn mice at the 3 hours and/ 
or 24 hours time points compared to the U3 sequences 
from no burn mice (Table 2 and Additional file 2: Table 
S2). Among the nine TREs are BRNF (Brn POU domain 
factors), AIRE (Autoimmune regulatory element binding 
factor), and NRF1 (Nuclear respiratory factor 1), which 
have been reported to be linked to various stress responses 
[25,26]. Furthermore, it appears that three TREs occurred 
more often in the U3 sequences of B-cells than T-cells 
whereas the U3 sequences from T-cells had a higher 
frequency of 16 other TREs in comparison to B-cells 
(Table 2 and Additional file 2: Table S2). Among the 16 
TREs, MAZF (Myc associated zinc fingers), GATA (GATA 
binding factors) and MEF2 (myocyte-specific enhancer 
binding factor) have been reported to be involved in T-cell 
activation [27-29]. A total of 122 TREs, including the 28 
TREs described above, were identified from the 100 unique 
sequences, and the comprehensive TRE profiles were 
summarized in Additional file 2: Table S2 and Additional 
file 3: Table S3. 



It is probable that the TREs {e.g., BRNF, AIRE, NRF1), 
which were mapped on the U3 sequences almost exclu- 
sively derived from burn mice, are closely associated 
with the modulation of MuLV-ERV expression in con- 
junction with a transcription environment dynamically 
formulated by burn-elicited stressors. 

In silico cloning and characterization of putative 
MuLV-ERVs harboring the unique U3 sequences 

To investigate the biological properties of the putative 
MuLV-ERVs harboring the U3 sequences identified in 
this study, the C57BL/6J mouse genome was surveyed to 
map the proviral loci using the 100 unique U3 sequences 
as mining probes. From the in silico mining experiment, 
49 ERVs were found with two identical LTR sequences 
at both ends and two, MuERV-49 and MuERV-26, 
retained slightly different LTR sequences due to a few 
point mutations. Forty-nine of the 51 putative MuLV- 
ERVs identified from this survey, including 15 full-length 
proviruses retaining a set of intact coding potentials for 
gag pol and env polypeptides, have been reported in 
previous studies [14,30-33]. The two newly identified 
MuLV-ERVs are: one full-length provirus (8,759 nucleo- 
tides) with intact coding potentials and one defective 
provirus (6,846 nucleotides) (Table 3). An analysis of the 
primer binding site revealed that both MuLV-ERVs had 
the tRNA Gln priming sequence. All 51 putative MuLV- 
ERVs were mapped using the U3 probes derived from 
both no burn and burn mice, except for BM-a-2.11e 
[31], which was identified using a U3 probe unique for 
B-cells of burn mice. On the other hand, 30 putative 
MuLV-ERV loci were retrieved using U3 probes exclu- 
sively isolated from T-cells. Future studies focusing on 
the biological properties of these putative MuLV-ERVs 
would be of interest: 1) effects {e.g., cytotoxicity) of 
their replication on immune cells in conjunction with 
infection followed by random integration into the 
genome and 2) functions {e.g., regulation of inflamma- 
tory mediators) of their gene products {gag, pol, and 
env polypeptides). 

Identification of genes neighboring the MuLV-ERV loci 

The U3 promoter and the enhancer sequences embedded 
in MuLV-ERVs can affect the transcriptional activities of 
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identified in the U3 sequences derived from specific experimental 
groups: burn group (green), B-cell group (yellow), and T-cell group (blue). 

Additional file 3: Table S3. Abbreviations: transcription regulatory 
elements. 

Additional file 4: Table S4. Annotated genes neighboring 100 Kb 
upstream and downstream of the individual MuLV-ERV integration loci. 
Chr (Chromosome number), Ori (Orientation). Gray highlight indicates the 
annotated genes harboring an ERV in the intron. *MuLV-ERV newly 
identified in this study. Please refer to Table 3 for the MuLV-ERV names. 



the neighboring genes via the promotion/repression of 
transcription, and alternative splicing and polyadenylation 
[34]. In this study, we identified the genes located within 
100 Kb upstream and downstream of each MuLV-ERV in- 
tegration site. Of a total of 131 genes identified, nine genes 
had MuLV-ERVs integrated into their introns (Additional 
file 4: Table S4). It will be interesting to examine how the 
expression patterns (depending on lymphocyte type or 
burn-elicited stress) of these proviruses are correlated with 
the expression of the relevant genes, for example, the rela- 
tionship between BM-a-2.11e MuLV-ERV, which was 
identified using a U3 probe uniquely identified from the 
B-cells of burn mice, and Camk2b (calcium/calmodulin- 
dependent protein kinase Il-beta) locus. The Camk2b ex- 
pression is associated with T-cell activation and cytotoxic 
T-cell activity [35,36]. It is possible that the burn-elicited 
stressors may affect the promoter and/or enhancer activ- 
ities of certain MuLV-ERVs paralleling modulation of the 
expression of adjacent genes in a specific lymphocyte type. 

Conclusion 

It is expected that burn injury elicits a range of stress 
signals that are directly linked to pathologic changes in 
the immune system and other organs [37]. The majority 
of the underlying mechanisms regarding how these stres- 
sors contribute to post-burn pathogenesis still remain to 
be characterized. The findings from this study suggest that 
certain MuLV-ERVs are differentially modulated depend- 
ing on the lymphocyte type, location of lymphoid organ, 
and duration after injury in conjunction with the burn- 
elicited stressors. As a result of activation of certain 
MuLV-ERVs following burn injury, the gene products 
and/or newly produced virus particles of the activated 
ERVs may participate in specific burn pathogeneses, such 
as immune disorder and cytotoxicity, in a cell type-specific 
manner. In fact, our laboratory has reported that some 
MuLV-ERVs are activated in the liver of mice after burn 
injury which resulted in evident hepatic damage [33]. Fur- 
ther investigations are needed to elucidate the link be- 
tween lymphocyte type-specific post-burn MuLV-ERV 
activation and systemic inflammatory disorder as well as 
the other aspects of burn pathogeneses. 

Additional files 
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